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Recent experimental observations have reported that, rather than randomly distributed, nitrogen 
atoms may prefer to be located on one of the two sub-lattices of graphene. It has been suggested that 
such a preference may present a possible avenue to tailor the band gap of graphene whilst maintaining 
its excellent electronic transport properties. Among the proposed mechanisms to explain this effect 
is the suggestion that long-range inter-impurity interactions mediated by the conduction electrons 
of graphene may give rise to the asymmetry between sub-lattice occupations. Electron-mediated 
interactions are known to be prevalent not only between N atoms but also between impurities that 
are adsorbed to a specific location within the hexagonal structure, namely to the top of the carbon 
atom. Furthermore, this interaction is known to become more long-ranged as the dimensionality 
of the system is lowered. For that reason, in this paper we investigate whether a similar sub¬ 
lattice asymmetry appears in the case of metallic carbon nanotubes doped with Hydrogen adatoms. 
Our results indicate that similar sub-lattice asymmetries are observable and even more pronounced 
in small-diameter CNTs with a dilute concentration of adsorbates, diminishing with increasing 
diameter or impurity concentrations. 


I. INTRODUCTION 

Since the discovery of Carbon Nanotubes (CNTs) in 
1991 1 there has been a huge amount of interest in ex¬ 
ploiting their unique mechanical and electronic proper¬ 
ties, with a wide range of applications such as energy 
storagd^Hl , composited® and electronic^®. The ex¬ 
traordinary properties of CNTs come from their atomic 
structure, being constructed of carbon atoms bonded via 
sp 2 bonds and forming a tubular hexagonal lattice - a 
quasi ID structure. It is well known that the introduction 
of dopants to nanostructures can change their electronic 
properties, and current research with doping CNTs aims 
to tailor the electr onic structure through controlling the 
dopant positioning^ 11 . 

Graphene, another novel carbon-based material, has 
had a similar level of research interest. While it has the 
same hexagonal structure as CNTs, it is totally flat - in¬ 
deed it can be imagined that CNTs are simply rolled up 
sheets of graphene. One common feature is that both 
materials, due to their hexagonal lattice structure, can 
be modelled as having two inter-penetrating triangular 
sublattices, is illustrated in Fig. [l] . Recent experiments 
of nitrogen doping in chemical vapour deposition-grown 
graphene have shown that dopants have a preference for 
occupying the same sublattice, forming large domains of 
many dopants on only one sublattice^®21 Even post¬ 
synthesis doping by direct ion implantation followed by 
heat treatment has shown preference for same sublattice 
configuration^^, suggesting the phenomenon can occur 
in other scenarios. Indeed, observations of same sub¬ 
lattice preference have been reported for low concentra¬ 
tion Molybdenum impurities 17 and high-coverage hydro¬ 
genated graphene 18 . Nevertheless there are many situa¬ 
tions where no asymmetry effect is observed, and more 
work needs to be done to clarify the situations in which 
one would expect to find the sublattice asymmetry in 


doping and to what degree. 

Despite this uncertainty, these findings have been par¬ 
ticularly interesting to the graphene community for many 
reasons. Firstly, it has been shown that this asymmetry 
in doping can be exploited in order to control the elec¬ 
tronic properties of the material introducing the possi¬ 
bility to open a band gap whilst preserving the excel¬ 
lent transport properties of its pristine state® 22 . Other 
work has investigated the novel properties that arise 
from edges and grain boundaries in the sample 23 and 
its potential use in ferromagnetic response from doped 
graphend^. Similarly, research with metallic CNTs 
doped with nitrogen has shown a large difference in trans¬ 
port response between segregated and unsegregated dop¬ 
ing pattern^ 10 . The question of whether it is experimen¬ 
tally feasible to induce a ferromagnetic state in either 
doped CNTs 2 - 5 27 and graphend^® hydrogeiPH, which 
arises from sublattice segregation of defects, can also be 
re-examined. 

The mechanism behind the segregation effect is de¬ 
bated and several theories have been proposed, for exam¬ 
ple the effect could occu r through nucleation during the 
CVD growth procesd^ 32 or alternatively from impurity 
interactions mediated by electrons in the host? 33 34 where 
it is known these interactions can be long-ranged 35 like 
the related Ruderman-Kittel-Kasuya-Yosida (RKKY) in¬ 
teraction for magnetic impurities 36 40 . The latter expla¬ 
nation captures the experimental observation of nitrogen 
asymmetry in both CVD synthesis and ion implantation 
with subsequent heat treatment, however the two mech¬ 
anisms may be complimentary to each other. 

We propose that if the second theory is correct that a 
similar effect should not only be observable in CNTs, 
but be more robust. Hypothetically the effect could 
occur in many metallic nanostructures but the mathe¬ 
matical formalism of the interaction and its similarity 
to the RKKY interaction suggests the effect should be 
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more pronounced in CNTs. The RKKY interaction is 
an effect by which two magnetic impurities can have 
their moments aligned or anti-aligned depending on their 
separatioiF^®, and this effect is mediated by conduc¬ 
tion electrons in the host. Additionally, the decay of the 
RKKY between two impurities separated by a distance 
(D) in any system follows its dimensionality (d), gener¬ 
ally going as the function D ~ d , hence in CNTs the decay 
is D~ x compared to a quicker D~ 2 in graphend^HIl. As 
both the RKKY and the proposed asymmetry effect arise 
from symmetry breaking in the lattice causing s o-call ed 
Friedel oscillations in the local density of stated^® it 
seems that CNTs are a natural candidate to test the lat¬ 
ter effect. 

In this work we will demonstrate that sublattice asym¬ 
metry of adsorbed atomic hydrogen atoms in CNTs 
should arise in a similar way to substitutional nitrogen 
dopants in graphene. Monoatomic hydrogen is perhaps 
the ’simplest’ dopant to consider from both an experi¬ 
mental and a theoretical point of view, and the doping 
of CNTs and other nanostructures with hydrogen is well 
understood due to the potential use for future energy 
storagd^EBHSI] Although nitrogen doped CNTs are well 
documented in the literature, hydrogen is an adsorbed 
impurity, instead of being an sp2 bonded substitutional 
impurity like nitrogen. The advantages of using hydrogen 
instead of nitrogen is that the doping can be applied post¬ 
synthesis of the graphene sheet and that the adsorbates 
can migrate atop the graphene with little energ}l^JMil 
This is not the case for nitrogeii 14 ^ 62 ! ^, although there 
is limited evidence to suggest that a subtle version of the 
sublattice asymmetry effect can be produced using high 
temperature annealing after post synthesis doping 1 ^. 

The mathematical framework developed in recent work 
for the asymmetry effect in graphene^ will be used as a 
basis for Monte-Carlo simulations where the additional 
effects of the geometry and size of the CNTs can be in¬ 
vestigated, along with the propensity for the adsorbed 
hydrogen to cluster. 


II. METHODS 

A. Tight-binding model and inter-impurity 
interactions 

While the mathematical framework for the as ymme try 
effect has been covered in-depth in recent work 3 ^^, an 
overview is presented here for self-completeness. Using 
a nearest-neighbour tight-binding method for graphene 
yields a Hamiltonian of the form H = — J2ij Uj describ¬ 
ing the hopping between nearest neighbour atoms i and 
j with the energy tij = 2.7eU. From this one can derive 
the Green Functions (GFs) for CNTs with both armchair 
(ACNT) and zigzag (ZZNT) geometries through inver¬ 
sion of the graphene Hamiltonian and enforcing certain 
periodic boundary conditions, a process explained in de¬ 
tail in the literature 4 -^^. As CNTs can vary in circum- 



FIG. 1. Schematic representation of a nanotube wrapped 
in the zigzag direction (ZZNT) doped with atomic hydrogen 
impurities, (a) shows three adsorbates bonded to different 
carbon host sites and colour coded depending on the sublat¬ 
tice of the host, also shown are the lattice vectors ai and a 2 . 
This is a close-up of the dashed box region in the larger ZZNT 
system shown in (c), where an ensemble of A — 15 impurities 
are shown in a central region of length L. Throughout the 
text we will refer to the degree of sublattice asymmetry in 
such a system using the percentage of impurities on the ma¬ 
jority sublattice. It is evident in this example case that there 
are 8 impurities on the white sublattice and 7 on the black 
sublattice, so the measured asymmetry is DE = 53.33% and 
the doping appears random. The central region is replicated 
along the axial direction, indicated by the ghosted copies to 
the left and right, in order to minimise finite-size effects. To 
exaggerate the concept of the central and replicated cells, each 
one is shown as a separate cylinder however in the model no 
such physical separation exists. The top-right schematic (b) 
shows the possible locations within the dashed blue circle re¬ 
gion that an impurity at site A in the central unit cell can 
move to if randomly selected in the iterative Monte-Carlo pro¬ 
cedure described in the main text, colour coded for opposite 
(green) and same (red) sublattice sites. 


ferential size, we follow the convention of ACNT n and 
ZZNT m to refer to armchair and zigzag tubes of circum¬ 
ferences defined by the lattice vectors ai and a 2 in Fig. 
|T] (a). If the circumference is defined as C = rqai + r 2 a 2 
it follows that a ANCT m (ZZNT m ) has a circumference 
defined by rq = m, r 2 = m (rq = m, r 2 = — m). 

Tight-binding is one of the more computationally vi¬ 
able methods for quantum mechanical Monte-Carlo sim¬ 
ulations, which rely on a huge amount of random sam¬ 
ples of a system to arrive at a reliable estimate of the 
ground state, due to its simplicity. Additionally, it is 
well known that the electronic structures of graphene 
and CNTs are described particularly well by the tight- 
binding method^. In the tight-binding regime the ad¬ 
dition of impurities and their effect on a pristine system 
with GF g is a relatively straightforward process by use 
of the Dyson Equation and a suitable description of the 
impurities. Atomic hydrogen adsorbates can be char¬ 
acterised thus by an onsite energy e a = 0.66£ relative 
to the Dirac point energy of graphene and higher than 
the system Fermi energy, and hopping integral r = 2.2 1 
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with a single carbon host site 67 . In the case of nitro¬ 
gen it was found that the exact parameterisation of the 
impurity can affect the quantitative results but that the 
qualitative behaviour remains the same. The parametri- 
sation of hydrogen used here captures the close range 
behaviour expected of pairs of hydrogens in graphene 
and CNTs whereby th ey prefer t o occupy opposite sub¬ 
lattices to each othei^^ 9 F 6Q * 68 * 69 [ The aim of this work 
is to model whether sublattice asymmetry can occur in 
nanotubes with adsorbed hydrogen impurities, so the be¬ 
haviour produced by this parameterisation should be suf¬ 
ficient for this purpose. Whilst numerous calculations of 
pairs of hydrogen adatoms, based on density functional 
theory (DFT), have shown that opposite sublattice con¬ 
figurations are more preferable, these are inadequate for 
considering the long-range disorder to be modelled here. 
Firstly, the symmetry breaking caused by the adsorb- 
tion of the dopant and the subsequent Ep shift will in¬ 
duce Friedel oscillations in the electron density leading 
to a long-range spatial oscillation in the energetics of the 
preferred configuration, and that it is the long-range be¬ 
haviour that leads to the asymmetry, as we will show. 
Also, the periodic nature of the DFT unit cell used to 
calculate these preferred configurations, effectively form¬ 
ing a superlattice, is also problematic when considering 
these interactions especially when a small cell is used. 

To calculate the most energetically favourable configu¬ 
ration for two or more impurities, the total energy change 
of the system must be calculated. When impurities are 
added to a system this energy can be calculated through 
the Lloyd Formula method^, and has the general form 

A E n = \j F dE Im Indet (i - g(E)V^j (1) 

where g is the GF of the pristine system, V is a matrix 
describing the perturbation introducing N impurities and 
Ep is the Fermi energy. 

For a system with two impurities the total system en¬ 
ergy change AE% can be separated into a relative separa¬ 
tion and configuration dependent term SE arising from 
the interaction of the impurities, and a term AEi which 
is the energetic cost of adding a single impurity in the 
system and is therefore position independent. As a re¬ 
sult it can be shown that 

SE = A E 2 - 2AEi (2) 

which will henceforth be referred to as the ” interaction 
energy” and has a long-range decaying oscillatory func¬ 
tional fornP^AE ~ cos 2 q f d w p ere Q F associated with 
the Fermi wavevector. The calculation of the ground 
state for a system with only 2 impurities is relatively 
simple computationally, but the complexity of this cal¬ 
culation increases with N\ for N impurities making the 
exact calculation prohibitively expensive. It is possible, 
however, to approximate the ground state using Monte- 
Carlo methods, and we demonstrate our methods for do¬ 
ing this in the following section. 


B. Monte-Carlo modelling of a finite concentration 
of impurities 

The strength of Monte-Carlo methods are that they 
can be used to simulate large complex systems where ar¬ 
riving at an exact numerical solution is difficult due to 
the scale or complexity of the problem. These techniques 
have previously been used to look at sublattice segrega¬ 
tion in general adsorbates on graphenP} and they will 
be used in this work as outlined below. 

Consider a section of nanotube of length L with a dis¬ 
tribution of N adsorbates randomly attached to host car¬ 
bon atoms in the system such as that in Fig. [l] Using 
the expression for SE and assuming only pairwise inter¬ 
actions, as higher order interactions will decay quicker, 
one can approximate the ground state by applying many 
iterations of choosing an impurity, summing all pairwise 
interaction energies between itself and the other impuri¬ 
ties, and comparing this energy to the same energy when 
the impurity has been moved to a nearby site (hereafter 
referred to as ^switch)? as demonstrated in Fig. [l] (b). If 
the second energy is lower then we assume that this low¬ 
ers the total system energy, and the impurity is moved. 
If the first energy is lower then no change is made. This 
technique is similar to that used for the computational 
approach to the Ising model of ferromagnetism in a 2D 
lattice. The effect of a Boltzmann Temperature Tp can 
also be introduced, whereby there is a finite probabil¬ 
ity ^ e ^^switch /k B T B m0 ve the impurity even though 
the first energy may be lower, where kp is the Boltzmann 
constant. It should be clarified that this iterative method 
of moving adsorbates is not the same as a dynamics cal¬ 
culation. While energy barriers for diffusion have been 
calculated to be on the order of < leV^HH] we choose to 
ignore these as the aim of this work is to find the ground 
state of the system and not model the physical diffusion 
of adatoms. Moreover, recent research has shown it is 
possible to make diffusion easier by doping CNTs with 
substitutional nitrogerP^. 

If a concentration of impurities (p) is added to the sys¬ 
tem it is natural to assume the number of electrons and 
therefore the Fermi Energy Ep will change. Ep is ad¬ 
justed with concentration according to DFT results for 
different length CNTs doped with a single hydrogen. It 
is found that the relationship is linear and approximated 
well up to p ~ 10% by assuming that there is ^ electron 
transfer to the host system per hydrogen and that the 
density of states profile is relatively unchanged. Other 
methods of finding this Ep shift have found a higher 
charge transfer from adatom to the graphene 7 ^ which 
would result in shorter oscillation periods and a reduction 
of the predicted segregation effect. 

To minimize finite-size effects which will be intrinsic 
to our methodology, it is necessary to impose some extra 
conditions. Firstly the central region and the impuri¬ 
ties are replicated along the axis of the nanotube once 
in each direction in order to approximate the effect of 
an infinite system. A schematic of this is shown in Fig. 
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[l] (c). When iterating, if an impurity in the central re¬ 
gion is moved then so are their ’twins’ in the replicated 
sections. Because of the oscillatory nature of the inter¬ 
action we impose a cut-off of one oscillation period of 
SE , beyond which the interaction is assumed to be zero. 
This is to ensure an even contribution to the interaction 
energy for all the impurities regardless of their position 
in the central region. Due to the scaling between length 
L, impurity concentration p and Fermi Energy Ep, this 
cut-off length is approximately the same as L and for 
the example system presented in Fig. [2] is approximately 
D = 180a, where a = 0.14nm is the carbon-carbon bond 
length of graphene, and this can be seen more clearly 
in Fig. [2] (a). Using the mathematical techniques de¬ 
veloped in previous work concerning Friedel oscillations 
and the sublattice asymmetry effect in nitrogen-doped 
graphene 33 34 1 it can be shown that the oscillation period 
changes as - p . 

The following calculations were done using approxi¬ 
mately N = 40 impurities for ACNTs and approximately 
N = 30 for ZZNTs, 107V 2 iterations and 250 randomly 
generated systems for each concentration and tempera¬ 
ture combination. The cut-off, directly linked to Ep and 
thus p, determines the exact value of L which in turn 
affects the value of TV, hence why there is a difference in 
N between ACNT and ZZNT systems. 


III. RESULTS 


A. Two hydrogen impurities 



FIG. 2. The interaction 5E between two hydrogen adatoms, 
with one fixed to a host at the origin and the second freely 
moved around to other carbon sites in the system, adsorbed to 
a ACNT6 host system with Ef — 0.033t. Although only two 
impurities are embedded in the system the Ef value has been 
shifted to simulate a doped system. The specific value corre¬ 
sponds to an impurity concentration of approximately 1.5% in 
the nanotube. Oscillations of 5E with distance along the ax¬ 
ial direction are shown in (a) for same sublattice (black) and 
opposite sublattice (red) configurations. The contour plot (b) 
extends this cross-section to sites in the circumferential direc¬ 
tion and the results are shown for same sublattice (left) and 
opposite sublattice (right) configurations where the short- and 
long-range behaviour is apparent. For clarity in the illustra¬ 
tion, the nanotube has been projected flat and the lattice sites 
are referred to by their circumferential and axial positions rel¬ 
ative to the carbon site hosting the hydrogen impurity at the 
origin. More detailed plots of the short-range behaviour in the 
left-hand side of (b), i.e. for same sublattice configurations, 
are shown for the same concentration (c), corresponding to 
the dashed box in (b), and a lower concentration of 0.1% (d). 


The behaviour of 5E for two hydrogens in a ACNT6 
is shown and discussed in Fig. [2j the most important 
characteristics to note are the difference in short- and 
long-range behaviours. When the impurities are close 
it is seen that they will prefer to occupy opposite sub¬ 
lattice arrangements, a behaviour occurring in graphene 
which was noted earlier in this work. This short-range 
behaviour gives way to long range decaying oscillations, 
complete with circumferential symmetry, and a distinct 
period 3 behaviour, characteristic of interactions between 
impurities with zig-zag separation in graphene-like sys¬ 
tems, can be observed 73 . The region of short-range be¬ 
haviour reduces with both increasing concentration, as 
shown by comparing plots (c) and (d), and increasing 
CNT circumference. The overall short-range and long- 
range behaviours and the oscillation period are very sim¬ 
ilar in ZZNTs with the exception that the period 3 be¬ 
haviour is in the circumferential direction, with long 
smooth decays along the axial direction. 


B. Monte-Carlo simulation of a finite 
concentration of hydrogen 

1. Finite Boltzmann temperature (Tb > 0) 


Monte-Carlo calculations were performed as described 
in Sec. |IIB| to investigate the appearance of sublat¬ 


tice asymmetry in ACNTs and ZZNTs with a range of 
sizes. The dependence of the sublattice asymmetry on 
the Boltzmann Temperature for a range of small diameter 
CNTs with ACNT and ZZNT geometries with different 
concentrations of hydrogen are shown in Fig. [ 3 ] A segre¬ 
gation of 100% corresponds to all impurities being on one 
sublattice whilst 50% corresponds to an equal distribu¬ 
tion between the two sublattices. As one would expect, 
increasing Tb causes the impurities to be randomly or¬ 
dered and the system will approach a 50% segregation 
(i.e. unsegregated or symmetric) state. 

The slower decay of the ZZNTs compared to ACNTs to 
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Boltzmann Temperature (K) 


FIG. 3. Degree of sublattice segregation, i.e. percent of 
impurities on one sublattice, vs. the Boltzmann temperature 
for small diameter ACNTs and ZZNTs with atomic hydrogen 
concentrations of 0.5% (black), 1.0% (red) and 1.5% (green). 
The relationship is analogous to that of magnetism and tem¬ 
perature in the Ising Model. It is apparent that increasing 
concentration leads to a lower expected segregation. This 
arises from the interplay between Ep and impurity separa¬ 
tion, both directly linked to the increased concentration of 
impurities, and this is discussed in detail in Sec. |IIIB2| 


the equilibrium state of 50% segregation originates from 
the difference in decay profiles and period 3 behaviour of 
SE. In ZZNTs it is harder to push an impurity out of 
an energetically stable configuration than in an ACNT, 
hence the broader temperature dependence for the former 
shown in Fig. [3j The relationship between segregation 
and Boltzmann Temperature indicates there is no meta¬ 
stable state at low temperatures, hence the dependence 
of segregation on the concentration can be investigated 
at Tb =0 and the result should be a good approximation 
of the ground state of these systems. 


2. Zero temperature (Tb — 0) 

The relationship between segregation and concentra¬ 
tion is plotted in Fig. [4] for a range of CNTs and the 
effect of increasing concentration and CNT diameter can 
be seen clearly, both leading to a reduction in the pre¬ 
dicted segregation effect. This can be understood as the 
long-range SE oscillations being the driving mechanism 
of the asymmetry, as was proposed for the effect for ni¬ 
trogen impurities in graphene 3 ^. Increasing concentra¬ 
tion or the CNT diameter reduces the average axial sep¬ 
aration of the impurities and the short-range behaviour, 
which favours opposite sublattice ordering, dominates. 
This can be corroborated by comparing the k-th nearest 
neighbour distributions of these systems and this is plot¬ 
ted in Fig. [5] comparing low (0.5%) and high (5.0%) im- 
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FIG. 4. Segregation vs concentration at Tb = 0 for CNT6 
(black), CNTg(red) and CNT 12 (green) of armchair (top) and 
zigzag (bottom) geometry, Monte-Carlo results are repre¬ 
sented by the dots and a line of best fit satisfying ~ 50 + p~ a 
for the fraction a ~ | is shown as a line. Increasing concen¬ 
tration or tube diameter leads to a reduction of the predicted 
segregation, tending to 50% in the limit of high concentration, 
a product of the reduction in the axial separation of impuri¬ 
ties leading to a dominance in the short-range behaviour of 
SE. 


purity concentrations, with segregated and unsegregated 
configurations respectively, and expected values for a ran¬ 
dom configuration in a ACNT6- Fig. [4] also suggests the 
asymmetry is generally expected to be higher in ACNTs 
than ZZNTs, where levels as high as 90% segregation 
should be observable at concentrations around 0.25% in 
a ACNT 6 . 

Following the trend of increasing diameter and decreas¬ 
ing segregation it seems a natural conclusion to expand 
the study to graphene by increasing the nanotube width 
to infinity, and where naively one should expect no seg¬ 
regation to be observed. This conclusion should not au¬ 
tomatically be made as there are subtle but important 
distinctions to be made between the case of nanotubes 
and graphene. Firstly, pristine graphene is semi-metallic 
but upon being doped by a finite concentration of hydro¬ 
gen the band structure is changed considerably, especially 
around the Dirac point ^1, leading to a more complex Ep 
shift and effects on the LDOS oscillations which drive 
the segregation effect. Secondly, these LDOS oscillations 
behave much differently in graphene decaying as D ~ 2 
and the directional isotropy of the coupling SE would al¬ 
low for more complete analytic calculations as has been 
shown for substitutional nitrogen impuri ties in graphene 
where a segregation effect was also found^^ 3 | 2Q t 34 l. 

Comparison of the nearest neighbour data in Fig. [5] 
demonstrates that clustering occurs in both high and low 
segregation systems, reflecting previous findings for hy¬ 
drogen pairs, both in Sec. |III A| and in the literature, that 
when in close proximity they prefer to occupy host sites 
on opposite sublattice J 6Q|68 i It has also been shown by 
Hornaeker et al. that a similar clustering effect can be 
observed in graphene 74 . At lower concentrations (Fig. 
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FIG. 5. Histogram data for 1st (top), 2nd (middle) and 
3rd (bottom) nearest neighbour separations for H pairs in a 
ACNT6 at 0.5% (left) and 5.0% (right) concentrations, com¬ 
paring measured data (black) to what one would expect for a 
random distribution (red). The discretised lattice makes iden¬ 
tifying the relative sublattice arrangement of the large peaks 
straightforward. Both low and high concentrations show a 
large peak corresponding to the nearest neighbour opposite 
sublattice arrangement. The data for low concentrations also 
shows a number of peaks at distances beyond a few lattice pa¬ 
rameters, which correspond to same sublattice arrangements 
in much greater amounts than the random data. At high con¬ 
centrations there are numerous peaks for close separations at 
much higher frequencies than the random data which points 
to clustering of the impurities. The 2nd nearest neighbour 
data shows two close peaks corresponding to opposite (left) 
and same (right) sublattice arrangements, and these also oc¬ 
cur in the 1st and 3rd nearest neighbour data. 


[5] - left) the most frequently found configurations are 
same sublattice ones, despite the short-range effect, due 
to the large average axial separation coupled with the 
long-range oscillatory nature of SE. This is evident by 
the location of the peaks in the low concentration data, 
occurring almost entirely for separations corresponding 
to same sublattice configurations. Additionally only a 
small deviation in the expected nearest neighbour profile 
(with the exception of the short-range peaks) is observed, 
along with the calculated mean and median values. 

These characteristics are in sharp contrast to the high 
concentration data, shown in Fig. [ 5 ] (right). In this case, 
there is a much larger degree of clustering, indeed almost 
all nearest-neighbours will be on the opposite sublattice, 
and there is a huge difference in the histogram profiles 
compared to the randomly calculated ones. The over¬ 
whelming frequency with which one finds opposite sub¬ 
lattice configurations at such short range lends credence 
to the long-range versus short-range/axial separation hy¬ 
pothesis, where the behaviour of SE acts so to cluster the 
impurities in opposite sublattice configurations at short 
range. 

A clearer illustration of the clustering effect and the 


short-range behaviour is shown in Fig. [6] where the prob¬ 
ability of finding impurities in certain short-range con¬ 
figurations is shown and compared for low and high con¬ 
centrations and the random data, aggregated from the 
nearest neighbour data of Fig. [5] For both high and 
low concentrations there is a large increase, 8% and 36% 
respectively, in the probability of finding a nearest neigh¬ 
bour same sublattice configuration over what one would 
expect in a random system. There is also a large increase 
in same sublattice next-nearest neighbour configurations 
- a consequence of the clustering of impurities. Most in¬ 
terestingly there appear to be preferential configurations, 
which one would expect from the profile of SE in Fig. [2j 
For example, there is a 4- or 5-fold increase in likelihood 
for the same sublattice configuration corresponding to 
the top-right most black site (position ~ 7.9a) and, due 
to lattice symmetry, equivalent sites such as the middle 
right-most black site at high and low concentrations. An¬ 
other preferential configuration is found for a separation 
of ~ 5.2a where the probability is approximately doubled. 
This is less pronounced than results seen in hydrogenated 
graphene^ where much higher concentrations and both 
sides of the monolayer were used, but is an interesting 
effect resulting from the inter-impurity interactions. 


IV. CONCLUSIONS 

In this paper a combination of tight binding Green 
Functions and Monte Carlo methods have been used 
to model sublattice segregation of hydrogen adsorbates 
in carbon nanotubes, motivated by recent findings with 
sublattice segregation and ordering of impurities in 
graphene. Understanding this effect is key to tailoring 
the electronic properties of graphene, nanotubes and per¬ 
haps other nanomaterials, along with potential to manu¬ 
facture of novel adsorbate ordered materials and inducing 
a magnetic response from graphitic nanosystems. 

It was shown in this work that the segregation effect is 
most pronounced for small diameter nanotubes with di¬ 
lute concentrations of impurities, and the effect decreases 
when either of these two parameters are increased. The 
inter-impurity interactions which mediate this interac¬ 
tion have a competing short- and long-range behaviour 
which tend to prefer opposite and same sublattice con¬ 
figurations respectively, and lead to certain preferential 
configurations and slight patterning in the lattice. In 
small diameter tubes with a dilute concentration of impu¬ 
rities the long-range behaviour along the axial direction 
dominates and the system prefers the segregated state. 
Increasing the concentration of impurities or the tube di¬ 
ameter reduces this axial separation and the short-range 
behaviour takes over, forcing the system into a state 
where the sublattices have approximately equal numbers 
of dopants. The short-range behaviour also leads to clus¬ 
tering and energetically favourable configurations of the 
impurities, a phenomenon which has been shown before 
both theoretically and experimentally. 
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(a) 


(b) 




FIG. 6. (a) Approximate probability for an impurity to 

occupy a host site, given an initial impurity indicated by the 
red square, for same sublattice (black) and opposite sublattice 
(red) arrangements at low (top number) and high (bottom 
number) concentrations in a ACNT6. Due to the symmetry of 
the system it is possible to allocate irreducible sectors, shown 
by the green shaded areas, which contain the information on 
all possible positions in the nanotube in order to simplify the 
data in this figure. To give some context for these findings, 
part (b) shows the stark difference between the probabilities 
shown in (a) and values that would be expected for a random, 
unrelaxed system. 


The current state of research with respect to gen¬ 
eral hydrogen doping of CNTj^ ^ 50 * 55 ! and other ID 
nanomaterialsEf^ an interest piqued due to their use 
for hydrogen storage for fuel cells, along with the capa¬ 
bility to di rectl y identify dopants and their respective 
sublattice d 62164175 1 , suggests that testing the findings of 
this paper and previous work in graphen^^, that inter¬ 
impurity interactions lead to a segregation effect, is well 
within current experimental limits. 
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